Stationary Nonlinear Schrodinger Equation on Simplest Graphs: 
Boundary conditions and exact solutions 
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We treat the stationary (cubic) nonlinear Schrodinger equation (NSLE) on simplest graphs. For- 
mulation of the problem and exact analytical solutions of NLSE are presented for star graphs 
consisting of three bonds. It is shown that the method can be extended for the case of arbitrary 
number of bonds of star graphs and for other simplest topologies such as tree and loop graphs. The 
case of repulsive and attractive nonlinearities are treated separately. 
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I. INTRODUCTION. 



The nonlinear Schrodinger equation has attracted much attention since from its pioneering studies in early seventies 
of the last century [l|-[3j. Such attention was caused by the possibility for obtaining soliton solution of NLSE and 
its numerous applications in different branches of physics. The early applications of NLSE and other nonlinear PDEs 
having soliton solutions were mainly focussed in optics, acoustics, particle physics, hydrodynamics and biophysics. 
However, special attention NLSE and its soliton solutions have attracted because of the recent progress made in the 
physics and Bose-Einstein condensates(BEC). Namely, due to the fact that the dynamics of BEG is governed by 
r^ ' Gross-Pitaevski equation which is NLSE with cubic nonlinearity, finding the soliton solution of NLSE with different 
'~^| confining potentials and boundary conditions is of importance for this area of physics. 

fT^ , Many aspects of soliton solution of NLSE have been treated during the past decade in the context of fiber optics, 

K*" photonic crystals, acoustics and BEG (see books [J]-[8J and references therein). Both, stationary and time-dependent 
Z^ NLSE were extensively studied for different trapping potentials in the context of BEG. In particular, the stationary 
PsJ ■ NLSE was studied for box boundary conditions 3, 'l(j| and the square well potential ll]-[l4]. 

In this paper we treat the stationary NLSE on graphs. Graphs are the systems consisting of bonds which are 



C ' connected at the vertices 15 1. The bonds are connected according to a rule that is called topology of a graph. 
^"^ ■ Topology of a graph is given in terms of so-called adjacency matrix (or connectivity matrix) which can be written as 



iOpolog 



. ^ ^ I 1 if i and j are connected, 

^. C,,=C,,^l . ,, . z,j = l,2,...,F. 

rN otherwise, 
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The linear Schrodinger equation on graphs has been topic of extensive research recently (e.g., see review 
and references therein) . In this case the eigenvalue problem is given in terms of the boundary conditions providing 



continuity and current conservation jl6| - 20 1 



Despite the progress made in the study of linear Schrodinger equation on graphs, corresponding nonlinear problem, 
i.e., NLSE on graphs is still remaining as less-studied problem. This is mainly caused the difficulties that appear in 
the case of NLSE on graphs, especially, for the time-dependent problem. In particular, the problem becomes rather 
nontrivial and it is not so easy to derive conservation laws 13311 . It should be noted that during the last couple of 



years there were some attempts to treat time-dependent 33l. |34| and the stationary 35|, |37[ NLSE on graphs. Soliton 



solutions and connection formulae are derived for simple graphs in the Ref. 33[. The problem of fast solitons on 



star graphs is treated in the Ref. 3J|. In particular, the estimates for the transmission and reflection coefficients 
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are obtained in the limit of very high velocities. The problem of soliton transmission and reflection is studied in 
[35|through the numerical solution of the stationary NLSE on graphs. 

Dispersion relations for linear and nonlinear NLSE on networks are discussed in |36| . More recent treatment of the 
stationary NLSE in the context of scattering from nonlinear networks can be found in the Ref . [37] ■ In particular, the 
authors discuss transmission through a complex network of nonlinear one-dimensional leads and found the existence 
of the high number of sharp resonances dominating in the scattering process. The stationary NLSE with power 
focusing nonlinearity on star graphs was studied in very recent paper [38J , where existence of the nonlinear stationary 



states are shown for (5— type boundary conditions. In particular, the authors of consider [38[ a star graph with N 
semi-infinite bonds, for which they obtain the exact solutions for the boundary conditions with a =/= 0. The properties 
of the ground state wave function is also studied by considering separately the cases of odd and even A^. We note 
that our work treats NLSE on simplest graphs with finite- length bonds, aiming at obtaining its exact solutions for 
some types of the boundary conditions. 

Despite the fact that considerable interest to NLSE and soliton transport on networks can be observed during last 
2 years, the problem is still very far from being studied comprehensively. In particular, detailed treatment of the 
boundary conditions and exact solutions for simplest topologies are not presented in yet in the literature. In this work 
we present exact solutions of the stationary NLSE for three types of simplest topologies, such as star, tree and loop 
graphs. 

Motivation for the study of NLSE on graphs comes from the different practically important applications such as 



soliton transport in optical waveguide networks [39|, soltion dynamics in DNA double helix [2l|-[23| and living systems 
[24J and other discrete structures | 25| . 
An important applications of NLSE on networks is Bose-Einstein condensation (BEC) and transport of BEC in 



networks. This issue that has been extensively discussed recently in the literature [26|-[32[. We note also that networks 
can be used a the traps for for BEC experiments. 

It is important to note that earlier the problem of soliton transport in discrete structures and networks was mainly 
studied within the discrete NLSE 25|. However, such an approach doesn't provide comprehensive treatment of the 
problem and one needs to use continuous NLSE on graphs. The aim of this work is the formulation and solution of 
stationary NLSE on simplest graphs such as star, tree and loop which can be considered as exactly solvable topologies. 

This paper is organized as follows. In the next section we will present formulation of the problem for the primary 
star graph consisting of three bonds. Section III presents derivation of the exact solution for stationary NLSE on 
primary star graph by considering repulsive and attractive nonlinearities. In section IV we discuss the extension of 
these results to the case of other simplest topologies such as tree graph, loop graph and their combinations. Finally, 
section V presents some concluding remarks. 

II. TIME-INDEPENDENT NLSE ON STAR GRAPHS 

The problem we want to treat is the stationary (time-independent) NLSE with cubic nonlinearity on the primary 
star graph. The star graph is a three or more bonds connected to one vertex (branching point). The primary star 
graph consisting of three bonds, bi, 62, 63, is plotted in Fig. 1. The coordinate, xi on the bond, 5i is defined from 
to Li, while for the bonds bk,k — 2,3 the coordinates, Xk, are defined from Li to Lk- At the branching point we 
have Xk = ii. In the following we will use notation x instead of a;^, (A; = 1, 2, 3). Then the time-independent NLSE 
can be written for each bond as 

- ij'J + Pj\^j\^^j = X^^„ P, > 0, J = 1, 2, 3. (1) 

Eq.([T]) is a multi-component equation in which components are mixed through the boundary conditions and conser- 
vation laws. Most simplest conservation is the current conservation. For each bond the current is defined as 



ik{x) = 2Im 



^'^^'^^ dx 



b, 




FIG. 1: Primary star graph consisting of 3 bonds connected at a vertex. 



This current should disappear at the ends of the each bond that can be written in terms of the following conditions: 



■0feU=z 



0, ^^^ \,=L,= J^kMLk), fc = 2,3. 
ox 



The following local current conservation condition (KirchhofF law) should be valid at the branching point 

Ji(ii) =.h{Li) +J3{Li). 

""""lac 



More detail analysis of the boundary conditions on graphs can be found in the Refs. [19 
In this paper we will consider the following boundary conditions: 



MLi) = A2ML1) = A^il^siLi), 



i^^^i^)-iiiiM^)-i^^i^M^) 



i_ a 
3 



\x=Lt 



aipi{Li), A2A3 ^ 0. 



In special case of ^2 = ^3 = 1, these boundary conditions reproduce those considered in the Refs. |16lll7|. 

III. SOLUTION OF THE TIME-INDEPENDENT NLSE ON SIMPLEST GRAPHS 

Consider the NLSE on the star graph presented in Fig 1. with the following boundary conditions (A is real; 
A2 = ^/^Pl, A3 = ^WWi, a = ) 



V'i(a;)U=o = 0, V'2,3(a;)U=L2.3 = 0> 

v^V'i(ii) = VW2ML1) = VhMLi). 

where the wave function is normalized as follows: 

3 






'xWdx^l. 



(2) 
(3) 
(4) 



(5) 



3 = 1""^ 



In the following we will consider Dirichlet boundary conditions at the end- vertices of a graph, while at the branching 
points the boundary conditions given by Eqs. ([3|) and (Hj) are to be considered. We note that in the left-hand side 
of Eq. m the derivatives are taken with the plus sign, while in right-hand the signs of the derivatives are positive. 
Furthermore, to make our notations different than those of the Refs. 3J, l38|, the conditions given by Eqs. (|3|, j2 
and their generalizations will be called ^"— type conditions (the boundary conditions considered in the Refs. [3^ 
are denoted by 5— and 5'). 



Solution of the NLSE can be written in the form 

V'j(a;) = exp(nj(x))/,(x), j = 1, 2, 3, (6) 

where fj (x) is a real function obeying the equation 

-/;+/3,/j = A^/,. (7) 

Inserting Eq. ([6]) into Eq. ([TJwe have 

-n^f,+i?fj-fj' + (3,f!^)^fj- (8) 

Taking into account Ec^. Q and requiring that real and imaginary parts of Eq. should be zero (|H]) we get jj = const. 
The following relations can be obtained from Eq. ([3]) 

exp{iji)^/%fi{Li) ^ exp{iJl)^/]3^f2{Ll) = exp{i^'i)y^f^{Li). 

The latter is valid only under the conditions 

71 = 72 = 73 = 7: 

It is clear that the functions /i, /2, /s should obey Eq. (HJ - ([5]). 

Exact solutions of Eq.© for finite interval and periodic boundary conditions can be found in the Refs. [9|,[l0|. Here 
we consider this problem for the network boundary conditions given by Eq. ^. Partial solution of Eq. ([7]) satisfying 
these boundary conditions can be written as 

fi{x) = Bisn {aix\ki) , 
f2{x) = B2sn {a2{x - ^2)1^2) , 
fsix) = Bjsn {a3{x - £3)1^3) , 



where sn{ax\k) are the Jacobian elliptic functions. [4] 
The general solution can be written as 

fb{x) = Bbsn (abX + Sb\h) . 

Inserting the last equation into Eq. ([7]) and comparing the coefficients of similar terms we have 

'2 



where ab = ±1. 



Bb ^ ab\j ^abkb, \^^al{l + kl), (9) 

Pb 



Using the relations 40 1 



/ sri^ {a{x ~ c)\k) dx = —7^ / \\ — dn^ [aix — c)\k)\ dx — 

J a k^ Ja 



— {b -a) —E [am {a{b - c)) \k] H —E [am, {a{a - c)) |fc] 



am{u + 2K{k)\k) — tt ^ am{u\k)^ 



E{mr ± ip\k) = 2nE{k) ± E{ip\k), 

and taking into account Eqs. (U)-® and ©, wc obtain the following nonlinear algebraic system with respect to aj 
and kj, j = 1, 2, 3: 

y^-Bisn (aiLilfci) = ^MB2sn (a2(ii - i2)|fc2) = V^B^sn (agl^i - ^3)1^:3) , (10) 



cn(o;iLi|fci)dn(aiLi|fci) r=^cn(a2(£i - ^2)1^) rf« ("2(^1 - -^2)!^) 

VPi vP2 

Bsas 



/^ 



en (a3(ii - ^3)1^:3) dn (asiLi - ^3)1^3) = 0, (11) 



R2 r2 R 



2-^1 + Ti-(^2 - ii) + Tr(^3 - ii) 



o2 n2 n2 

= 1 + jj^E [am(aiii|fci)] + -^^ [am(a2(£2 - ^1)1^2)] + jj^E [am(a2(i3 - ^Ol^s)] • (12) 

Kfai K2Q;2 KgOfs 

Here i4r(fc) and E{k) are the complete elliptic integrals of first and second kind, respectively. 

In general case this system can be solved using the different (e.g., Newton's or Krylov's method) iteration schemes. 
However, below we will show solvability of this system for two special cases. 
First case. Let ^^ = ZT-I] ^ ^^^^re 711,2,3 e N U {0}. 
Choosing 

(Xl = K[ki), Q2,3 = 7 J-K(k2,3), 

we have 

— ai ~ a2 = as = a, fci = ^2 == fcs = k, (Ji ~ (72 — 0^3 ^ 1- 
It is clear that under these conditions Eqs. pUj) and ([TT]) are valid. Furthermore, it follows from Eq. p^ that 

9{k) ^ 2^^^^ (|i + ^^ + ^^) Kik) (Eik) - (1 - k^)Kik)) -1 = 0. (13) 



The solvability of Eq. ([131) is equavalent to that of NLSE on graph. Therefore we will prove the solvability of this 
equation. Since the following relations are valid 

lini g{k) = —1, lini g{k) = +00, 

and g{k) is a continuous function of k on the interval (0; 1) it follows that Ea. (|T3)) has a root. 

Second case. Now we show that there exist another solution of Eg. p^ . For the case when ai — 
i-ir^P+2n,Kik,) ^ a2,3 = (~i)"''-^^P+2^2.3^(fc^.3) ^ ^hcrc < p < 2i^(fci,2,3), 7ii,2,3 € N 7^1,772,713 cannot be odd 
or even at the same time, from Eqs. ^ and ([TU| we can obtain 

ai = a2 = «3 = ce, ki = k2 = kz = k, ui — 1, 02= dz — — 1- 
From Eq. (jlip we have 

(_l)ni (_l)n. (_i)n3 _ 

/3l /32 /33 



Furthermore, it follows from the last equation and Ea. p2|) that 



g{k) = 4 



(-l)"ip + 2niX(fci) 



Li 



f + ^ + ^]iKik)-Eik))-1^0. 

Pi P2 P3 , 



(14) 



We have 



lim g{k) ~ —1, lim g{k) = +oo. 



Since g{k) is a continuous function of k on the interval (0; 1), it follows from the last relations that Eg. ((H)) has a 
root. 

A. Case of attractive nonlinearity 

NLSE with attractive nonlinearity is of importance for a number of problems of BEC physics llOll. Earlier, it was 
treated for the box-boundary conditions in the Ref . 10 1 . Here we extend the results of the Ref . [lO| to the case of 
the star graphs. Consider also the case of NLSE with attractive nonlinearity given as 



^y - /3,#,|V, = A^V',, /3, > 0, J = 1, 2, 3. 



(15) 



We write the solution of Eq. (jlSp as 



ipj (x) = Bjcn {ajX + 5j\kj) , 
where the parameter 5j can be determined from the boundary condition given by Eq. ^ as 

5i = (2mi + l)K{ki),5j = (2mj + l)K{kj) - a^Lj, j = 2, 3, 
where mi, 1112,1113 E Z. The coefficient Bj obeys Eq. ([9]) with 

X^ = {l-2k^)a^j. (16) 

Taking into account Eqs. (ITU | - P^ and (|16p we obtain the following nonlinear algebraic system with respect to 



Xj, 



kj J = 1,2,3: 



/p^Bicn (aiLi + (2mi + l)X(fci)|fci) = y%B2cn (03(^1 - L2) + (2m2 + I)i^(fc2)|fc2) = 

V^Bgcn (a3(ii - L3) + (2to3 + l)K{k3)\k3) , (17) 

(1 - 2kl)al = (1 - 2kl)al = (1 - 2kl)al, (18) 



Bia 



Isn {aiLi + (2toi + l)K{ki)\ki) dn [aiLi + (2mi + l)K{ki)\ki) - 



i?2a2 

7^ 



V3 

sn {a2{Li - L2) + (2m2 + I)is:(fc2)|fc2) dn (a2(ii - ^2) + (2m2 + I)i^(fc2)|fc2) - 






sn (a3(ii - ia) + (2m3 + l)X(fc3)|fc3) dn (03(^1 - L3) + (2m3 + l)if (^3)1^3) = 0, 



(19) 



aik\ 



[E {am{aiLi + (2mi + l)K{ki)\ki)) ~ (2mi + l)£'(fci)] 



a2fc| 



[(2m2 + 1)-E(A:2) - E {am{a2{Li - L2) + (2m2 + l)K{k2)\k2))] 



[(2m3 + l)E{k3) - £; (am(a3(ii - L3) + (2m3 + l)X(fc3)|fc3))] 



= i + Bi 



1-kl 



Li + B' 



1 k2 

•^2 



(L2 - il) + S3" 



2 1 ^3 



Kg 



(is-ii). 



(20) 



Now we again consider two special cases . 

First case. Let ^^iL+i ^ f^^^^^, where toi,2,3 e N U {0}. 

Then choosing 

2mi + l 2m2,3 + l„., . ,„^. 

ai = r -f^(«i)> "2,3 = 7 J-K[k2,3)-, (21) 

^1 -(^2,3 — -'^1 

from Eqs. (|17p and(fT8|) we can obtain 

— ai = a2 = as = a, ki = k2 = k^ — k, cti = (T2 = fs = 1- (22) 



It is clear that under these conditions Eq. (1191) is valid. Furthermore, it follows from Eg. ((221) that 

gik) . 2i?!^^ (|i + ^^ + ^^) i^(fc) (i^(fc) - (1 - k^)Kik)) -1 = 0. (23) 

Again, the solvability of Eq. (P51) is equavalent to to that of NLSE with attractive nonlinearity. 
Therefore we will prove the existence of the roots of this equation. We have 

lini g{k) = —1, lim g{k) = +oo. 

fc— !-0 ' fe— S-l 

Then the continuity of function g{k) on the interval (0; 1) implies that Eq. (P5)) has a root. 

Second case. Now we show that there exist another solution of Eq. (f23| . Let ai = ^^ ni- mi- ) — UJ.^ q,^ ^ = 

_ TP+(4n..3-2^_.^3-l)i^(fc..3) ^ ^j^g^g < p < 2if (fcl,2,3), ^1,2,3 G N 4ni,2.3 > 2mi,2,3 + L 

Then we have 

|ai| = \a2\ = |a3| = a, ^1=^2 = ^:3 = k, ai = 1, (T2 = 0-3 = -1- (24) 

From the condition (fT9|) we have 

1 1 1 

Furthermore, from the last equation and normalization condition we have the following relation: 
±p+(4ni -2mi- l)K{k) 



g{k) = 2 



Li 



An,-2m,-l , 4n2-2m2-l , 4n3 - 2m3 - 1 , ^^^^^ _ ^^^^^ _ ^ ^ ^^ 



/3i /32 /33 

(25) 

Therefore, according to the asymptotical relations given by 

lim5(fc) = — 1, \iuY g(k) — +00 ^ (26) 

and from the continuity of function g[k) in the interval (0; 1), we can conclude that Eq. ()25p has a root. 

IV. OTHER SIMPLEST GRAPHS 

The above results can be extended for other types of graphs, such as general star graphs, tree (Fig. 2), loop (Fig. 3) 
and their combinations with appropriate boundary and vertex conditions. This can be done using the same approach 
as that in the Ref. [33]. 

A. Tree Graph 

Here we will give details of such extension for the tree graph. 




FIG. 2: Tree graph. 



In the following we assume that 5" ~ type boundary conditions are hold at the branching points, while at the end 
vertices the wave function becomes zero. 

Furthermore, we seek for the solution of Eq. ([7]) on the each bonds in the form 

fh = Bbsn (afcX + Sb\h) , 

where Sb are parameters that can be determined by the given boundary conditions: 6i = 0, Suj = —aujLuj. 

Again, as it was done in the previous section for star graph, one can show posibility of the exact solution for two 
special cases. 



The first case is given by the relations 

4ni + 1 _ 4(n^P - n'-^^) _ Anuj + 1 



Li Lii — Li 

(2) ^^(1) 



Jlij — IJli 



, 1 = 1,2, j = 1,2,?. 



where ni,nii,niij G N U {0}, nj^j-^ > n]^^' and 

Li Lii — Li Liij — Lii 

It is easy to show that the following relations are valid: 



Anuj + 1 






K{ku) 



ai = an = aiij = a, ki = ku = kiu = k, ai = an = 1, auj = -1. 



Furthermore, we have 



,(„,2(±!L±l)Mi. 



^+i: 



'-i 1ft .=, 



Lii — Li Y^ Liij — Li 



Pi 






K{k){K{k)~E{k))~l={). 



Solvability of the last equation is obvious. 
The second case is given by the relations 



ai 



-(-i)"ipi+2»i-ft:(fci) 



Li 



, OLU 






nij 



_ -(-l)"i'Jpii+2niij,g(fci.,) 



Su 



(-1)"^ piLii-(-l)"i. ptiLi+2in[\^ Lti-n['l'> Lt)K{kii) 
Lii-Li ' 



where < pi < 2K{klA^), < pu < 2K{ku,uj). 
From vertex conditions we obtain 



and 



ai = ait = aiij = a, ki = ku = ku^ = k 



(_i)m ^ (-1)"^-' _^^ (-1)""' ^ y^ (-ir- ^^_ 



A ^ A. 



^ /3iij 



It follows from the normalization condition that 



,^,,.4(tlr^l±^)[l 



- + E 



(2) (1) 3 



^S^+E 



/3; 



1 ^Hj 



{K{k)~E{k))- 1 = 0. 



Solvability of the last equation follows from the properties of the function g{k). 

B. Graph With Loops 

Similar results can be obtained in for the loop graph plotted in the Fig. [31 

The solution of Eq. (I?]) for this topology can be written in terms of the Jacobian elliptic functions as 

fi{x) = Bisn {aix\ki) , 

fj{x) = Bjsn {ajX + Sj\kj) , j = 2, 3, 

f4{x) = B^sn {a4{x — L4)|fc4) . 



fti 



^I^> 



FIG. 3: Loop graph. 



Again, two special cases will be considered. 
For the first case we put 






,(1),(2) 



l^^l^l^f , where 7^1,4, n^|^'^^^ G NU {0}, n^^^J^ > n^"^. Choosing 



.(2) ^ „(i) 



^t (2) (l)■^ 

"1^4 = ^1.4 K{ki^4), ^2,3 = L, a-Li ^(^2,3), 



<52,3 = 



4{4]lL2,3-ni%Li) 



+ 1 /V(fc2,3). 



L2.3 — L1 

We have aj = a, fcj = fc, j = 1, 2, 3, 4, cri = (72 = fs = 1, <T4 = — 1. From the normalization condition we get 
_ (4711 + 1)2 fii , L2-L1 L^-Li L4 



g{k) = 2- 



a 



Pi ii2 



The second case is given by 



, -^\K(k){K{k)-E(k))-\ = 0. 

P3 Pi 



,(2) ^(1)^ 



i-ir^Pi + 2n,K{ki) i-ir-''P2 - i-ir-^-Pi + 2(n^^^ - n'^i)K{k2,3) 
OL\ = , a2,3 = 



a4 = 



(-l)'>2 + 2n4if(fc4) 



-^2,3 ^ -^1 



^2.3 = 



(_l)"2,3p^L2^3 - {-ir-^'P2Li + 2(4^^L2^3 - n^2.iLl)K{k2^3) 



^2,3 



il 



10 

where < pi < 2K{ki,2,z), < p2 < 2K{k2.3,4)- 

Then from the 5"— type vertex conditions we obtain aj = a, kj — k,j — 1,2,3,4 and 

(-1)"^ (-iK' (-iK' ^n Lull , tilll , Li)!l = n 

/3i (32 133 ' (32 Ps 134 

Furthermore, according to the normahzation condition one can obtain the fohowing equation 

In both cases solvabihty of the equation g(k) — foUows from the asymptotical relations, liuik^Q g(k) — 
— 1, linifc^i g{k) = +cx) and from the continuity of the function g{k). 

V. CONCLUSIONS 

In this paper we have treated the stationary nonlinear Schrodinger equation with cubic nonlinearity for simplest 



graphs. Unlike the previous studies [33|, |3J, |37| , the lengths of the bonds are considered as finite. Therefore, our work 



can be considered as an extension of the results by L.D. Carr [9|, |lO|, [l^ to the case of networks. Formulation of the 
problem and its detailed treatment including method for finding exact solutions are presented for the case of star 
graph consisting of three bonds. The method is extended for other simplest topologies such as tree and loop graphs. 
The cases of attractive and repulsive nonlinearities are treated separately. The above results can be useful for the 
problems of BEC on networks and discrete traps, soliton transport optical waveguide networks, soliton excitation in 
DNA double helix, energy transfer in nanoscale networks etc. In principle, the method developed in this work can be 
extended for more complicated topologies, too. Currently such a study is on progress. 
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